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Abstract. Recently available near-infrared interferometric data on late-type stars show a strong increase of diameter for asymp- 
totic giant branch (AGB) stars between the K (2.0 - 2.4 /jm) and L (3.4 - 4.1 fjm) bands. Aiming at an explanation of these 
findings, we chose the objects a Orionis (Betelgeuse), SW Virginis, and R Leonis, which are of different spectral types and 
stages of evolution, and which are surrounded by circumstellar envelopes with different optical thicknesses. For these stars, we 
compared observations with spherically symmetric dust shell models. Photometric and 1 1 yum interferometric data were also 
taken into account to further constrain the models. - We find the following results. For all three AGB stars, the photosphere 
and dust shell model is consistent with the multi-wavelength photometric data. For a Orionis the model dust shell has a very 
small optical depth (0.0065 at 11 yum); the visibility data and model in K and L are essentially entirely photospheric with no 
significant contribution from the dust, and the visibility data at 1 1 /jm show a strong dust signature which agrees with the 
model. For SW Virginis the model dust shell has a small optical depth (0.045 at 1 1 pm); in K the visibility data and model are 
essentially purely photospheric, in L the visibility data demand a larger object than the photosphere plus dust model allows, 
and at 11 fim there was no data available. For R Leonis the model dust shell has a moderate optical depth (0.1 at 11 //m); in 
K and L the visibility data and model situation is similar to that of SW Vir, and at 1 1 /jm the visibility data and model are 
in agreement. - We conclude that AGB models comprising a photosphere and dust shell, although consistent with SED data 
and also interferometric data in K and at 1 1 yum, cannot explain the visibility data in L; an additional source of model opacity, 
possibly related to a gas component, is needed in L to be consistent with the visibility data. 

Key words. Techniques: interferometric - Radiative transfer - Circumstellar matter - Infrared: stars - Stars: late-type - Stars: 
individual: a Orionis, SW Virginis, R Leonis 



1. Introduction 

Supergiant stars and stars on the Asymptotic Giant Branch 
(AGB) are surroiinded by a dust envelope coupled with gas; 
for example, se e iHabinel ( 1 19961) for a review of these ob jects, 
IWillsonlll2000ll for a review of mass loss, and lTsuiil ( l2000i) for 
evidence of H2O shells around supergiants. 

AGB stars are subject to pulsation and other surface phe- 
nomena (shock waves, acoustic waves etc., see Lafon & 
Berruyer (1991)) that levitate matter above the stellar photo- 
sphere. Dust grains appear in the envelope at altitudes where 
the temperature falls below the condensation temperature. The 
radiation pressure drives the dust grain overflow far from the 



Send offprint requests to: Peter Schuller, e-mail: 
pschullerOcf a. harvard. edu 



star up to distances of a thousand stellar radii causing mass- 
loss. 

The mechanism for mass loss in supergiants is not so clear. 
The main difference compared to Mira stars is that due to 
the larger luminosity of supergiants the atmosphere is much 
more extended and has a much greater pressure scale height. 
Therefore, the mass loss may have a different driving force al- 
though pulsation may also play a role. 

Both kinds of stars are characterised by their large diameter 
and strong luminosity and are consequently very good candi- 
dates for infrared interferometry. In the next section, we briefly 
present the instruments which were used to acquire the interfer- 
ometric data discussed in this paper, namely FLUOR in the K 
band (2.0 - 2.4 yum) and TISIS in the L band (3.4-4. 1 yum). We 
selected three particular stars of three different types for which 
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we have a substantial set of data: the supergiant a Ori, the semi- 
regular SW Vir, and the Mira-type star R Leo. Furthermore, 
their envelopes have quite differing optical thicknesses. It ap- 
pears that for the objects with large optical depth (SW Vir and 
R Leo) the observed diameters in L are considerably larger than 
the respective K band diameters (see Sect. 12.31 1. Measurements 
in the L band are more sensitive to cooler material above the 
stellar photosphere than are the K band data. The question is: 
what is the physical nature of these layers? One possibility 
might be astronomical dust. Dust reprocesses the star radiation, 
in particular in the infrared. Depending on the dust abundance, 
the star's L band radiation might be superposed by a dust con- 
tribution that probes the inner region of the dust envelope at a 
few stellar radii away from the star. The measured L band ra- 
dius might consequently be increased compared to the K band, 
which is mainly fed by stellar emission. 

In Sect.|3] the modelling of the dust envelopes around these 
stars by a radiative transfer code is explained and the param- 
eters of the model are described. In Sect. |4] we compare the 
model output with the interferometric and photometric data. 
Finally, in Sect.|5]we discuss the limits of such models and the 
questions they raise about the L band observations. 

2. Observations and data reduction 

2.1. Instruments used 

The stars discussed here were observed with the IOTA 
(Infrared-Optical Telescope Array ) interferome ter operating 
(at that time) with two telescopes (Traub 1998). IOTA is lo- 
cated at the Smithsonian Institution's Whipple Observatory on 
Mount Hopkins, Arizona. It offers multiple baseline observa- 
tions providing visibilities at different spatial frequencies. This 
enables a model-fit to the visibility data to obtain information 
on the spatial intensity distribution of the sources. The beams 
were combined with FLUOR (Fiber Linked U nit for Optical 
Recombination) in the K band ( Coude du Foresto et aL 1998 ) 
and with its extension to the L band TISIS (Thermal Infrared 
Stellar Interferometric Set-up, Mennesson et al. 71999")). Beam 
combination with FLUOR is achieved by a single-mode fluo- 
ride glass triple coupler The fibres filter the wavefronts corru- 
gated by the atmospheric turbulence. The phase fluctuations are 
traded against photometric fluctuations which are monitored 
for each beam to correct for them a posteriori. The accuracy 
on visibility estimates measured by FLUOR is usually better 
than 1% for r nost sources (Perrin et al. 1998.) and can be as 
good as 0.2% llPerrinl2003h . 

TISIS is the extension of FLUOR to the thermal infrared 
(3.4 - 4.1 fim). A single coupler is used as beam combiner 
and photometric signals are not monitored. Approximate cor- 
rection of turbulence-induced flux fluctuations is achieved with 
the low-frequency part of the interferometric signal. Accuracy 
in the L band is therefore not as good as in the K band, but still 
better than without spatial filtering by fibres. The acquisition 
protocol is also different from that of the K band, since ground- 
based interferometric observations encounter a new difficulty at 
3.6 ;um. Thermal background emission produces a fluctuating 
photometric ofl'set which must be correctly subtracted in order 



Table 1. Some basic parameters of the studied obje cts taken 
from the General Catalogue of Variable Stars ( Kholo pov et alJ 
'1998) . The distances come from the HIPPARCOS catalogue 
(.Perrvman etalJl997i) . 



Star 


Spectral 


Variability 


Period 


Distance 




Type 


Type 


(days) 


(pc) 


a Ori 


Mil 


SRc 


2335 


131 


SW Vir 


M7III 


SRb 


150 


143 


RLeo 


MSIIIe 


Mira 


310 


101 



to achieve a good calibration. Hence each source observation 
must be bracketed by a sky background measurement (chop- 
ping) that reduces the instrument efficiency. 

2.2. Data selection 

We chose to include three stars in our study which span 
different spectral types and stages of evolution, namely the 
supergiant a Orionis (Betelgeuse), the semi-regular variable 
SW Virginis, and the Mira star R Leonis. Some basic param- 
eters of these objects are compiled in Tabled These stars are 
also some of the m ost observed in our sam ple. The L band data 
are extracted from lChagnon et aP ( l2002l) - we have used the 
data from the February-March 2000 observations. The K band 
data for a Orionis were presented in Perrin et al. (2003a) and 
observations cover the period February 1996 to March 1997. 
The K band data for SW Virginis were collected in May 2000 
and are presented in Perrin ( 2003 ). Lastly, the K band data for 
R Leonis were published in lPerrin et alJ (1 19991) and were col- 
lected in April 1996 and March 1997. Except for SW Virginis, 
there is a large time difference between the dates of the L band 
and K band data acquisitions. We believe this to not be a major 
issue - although these stars are known to be variable, the ef- 
fects we want to analyse are slower than the periodic changes 
of stars and larger in amplitude. 

In order to get the best constraints on the dust shell model 
spatial distribution, we have included interferome tric data ac- 
quired with ISI (Infrared Spatial I nterferometer , iHale et alJ 
ilOOd) ) at 1 1 . 1 5 yum as published bv lDanchi etaP (Il994l) . 

A journal of all used data points is given in Table|2l 

2.3. Uniform d/sk diameters 

The visibilities V,, and their standard deviations cr,, in the K and 
L bands were fitted by a uniform disk model by minimising the 
quantity: 



X 



=2 



(Vi-M(cf>iju,Si)f 



(1) 



/■=! 



where n is the number of available measurements in the respec- 
tive band. M is the uniform disk model: 



2Ji (jTipiJDS i) 



(2) 



with S i the spatial frequency and 0ud the uniform disk diame- 
ter. Ji is the Bessel function of first order. 



p. Schuller et al.: Dust shell models and interferometric data of late-type stars in the near-infrared 



3 



Table 2. Journal of interferometric data. K band (FLUOR) da t a for a Orionis have been presented in lPerrin e t alJ (l2003a(), for 
SWVirginis in|Perrin (2003), and for R Leonis in |PerrinetalJ ( ll999l) . All L band (TISIS) data were taken from Chasnon e t alJ 
("2002"). The ISI data points for a Orionis and R Leonis were extracted from lPanchi et alJ 11994) using the DEXTER online tool 
at .AD&c2003.) . No exact position angle was given along with the ISI data, therefore we indicated the basic geographic orientation 
of the baselines used. Where known for the objects, the phase CD of the visual variability was indicated. 



Object 


Spatial Frq. 


Visibility 


Pos. Angle 


Object 


Spatial Frq. 


Visibility 


Pos. Angle 




(10^ rad"') 




(deg) 




(10*^ rad"') 




(deg) 


a Ori 


FLUOR 






R Leo 


FLUOR, Apr. 1996, O = 0.24 






46.41 


0.2232 ± 0.0018 


85.91 




65.76 


0.2723 ± 0.0040 


86.9 




46.82 


0.2045 ± 0.0020 


80.47 




65.76 


0.2718 ± 0.0043 


92.5 




48.14 


0.1847 + 0.0024 


71.96 




66.38 


0.2723 ± 0.0041 


80.3 




50.39 


0.1253 +0.0016 


63.01 




69.82 


0.2373 ± 0.0033 


65.4 




88.30 


0.0949 ± 0.0047 


69.89 




71.57 


0.2127 ± 0.0035 


60.7 




90.92 


0.0843 ± 0.0024 






FLUOR, Mar. 1997, O = 0.28 






133.99 


0.0510 ±0.0010 


79.06 




47.94 


0.4556 ± 0.0041 


86.0 




TISIS, March 2000 






50.08 


0.4248 ± 0.0036 


68.5 




35.99 


0.517 ±0.015 


88.85 




66.89 


0.1934 ± 0.0036 


77.0 




35.99 


0.528 ± 0.014 


89.49 




68.27 


0.1751 ±0.0029 


70.9 




36.34 


0.521 ±0.011 


81.08 




70.69 


0.1656 ±0.0050 


63.2 




53.32 


0.085 ± 0.001 


108.69 




70.96 


0.1761 ±0.0071 


62.5 




53.71 


0.081 ±0.001 


110.40 




91.07 


0.1255 ±0.0096 


71.2 




ISI, Oct. 1988 


- Oct. 1989 






97.21 


0.1188 ± 0.0060 


58.5 




1.64 


0.54 ± 0.04 


90 




150.62 


0.0770 ± 0.0032 


63.1 




1.82 


0.57 ± 0.02 


90 




TISIS, Mar. 


2000, O = 0.81 






2.24 


0.61 ± 0.01 


90 




39.29 


0.515 ± 0.011 


168.48 




2.36 


0.65 ± 0.01 


90 




39.46 


0.528 ± 0.011 


67.44 




2.51 


0.62 ± 0.02 


90 




40.35 


0.504 ± 0.013 


62.78 




2.75 


0.59 ± 0.01 


90 




52.76 


0.249 ± 0.002 


86.66 




2.96 


0.59 ± 0.01 


90 




52.91 


0.253 ± 0.003 


96.07 




3.26 


0.56 ± 0.02 


90 




53.69 


0.258 ± 0.003 


103.86 




3.53 


0.60 ± 0.01 


90 




53.88 


0.245 ± 0.002 


105.20 




ISI, Oct. - Nov. 1992 






53.94 


0.251 ± 0.003 


105.55 




6.91 


0.58 ± 0.01 


123 




54.04 


0.246 ± 0.002 


106.15 




7.39 


0.59 ± 0.01 


123 




54.10 


0.248 ± 0.002 


106.58 




7.57 


0.56 ± 0.02 


123 




54.12 


0.227 ± 0.002 


106.74 




7.96 


0.54 ± 0.02 


123 




54.85 


0.229 ± 0.003 


110.85 




8.67 


0.56 ± 0.02 


123 




ISI, Oct. 1988, * ~ 0.45 






ISI, Nov. 1989 


- Oct. 1992 






2.44 


0.86 ± 0.05 


90 




8.97 


0.61 ± 0.02 


113? 




2.54 


0.83 ± 0.04 


90 




9.18 


0.58 ±0.01 


113 




ISL Oct. 1990,0 = 0.81 






9.30 


0.60 ± 0.02 


113 




7.30 


0.55 ± 0.06 


113 




9.78 


0.56 ± 0.03 


113 




7.51 


0.63 ± 0.03 


113 




9.96 


0.55 ± 0.02 


i 1 J 




8.02 


0.53 ± 0.04 


113 




10.22 


0.54 ± 0.02 


113 




8.23 


0.65 ± 0.05 


113 




10.55 


0.56 ± 0.02 


113 




8.62 


0.59 ± 0.03 


113 




10.70 


0.59 ± 0.02 


113 




9.10 


0.54 ± 0.03 


113 




10.91 


0.54 ± 0.01 


113 




9.30 


0.65 ± 0.06 


113 




10.97 


0.58 ± 0.02 


113 




9.61 


0.53 ± 0.03 


113 




11.30 


0.54 ± 0.01 


113 




9.81 


0.59 ± 0.03 


113 




11.45 


0.48 ± 0.02 


113 




10.02 


0.60 ± 0.02 


113 




11.57 


0.55 ± 0.02 


1 1 Q 
i 1 J 




10.20 


0.56 ± 0.03 


1 1 ^ 

Hi 




11.69 


0.50 ± 0.02 


113 




10.32 


0.57 


113 




11.81 


0.52 ± 0.01 


113 










SW Vir 


FLUOR, May 2000, O = 0.75 














64.29 


0.697 ± 0.010 


61.31 












81.89 


0.522 ± 0.014 


93.93 












81.93 


0.527 ± 0.014 


94.23 












TISIS, Mar. 2000, O = 0.35 














32.63 


0.830 ±0.012 


76.01 












32.73 


0.848 ±0.012 


75.39 












46.80 


0.706 ± 0.008 


98.57 
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Table 3. Uniform disk radii in the L and K bands. 



Star 








l/(Pk 




(mas) 


(mas) 






a Ori 


21.37 ±0.01 


21.63 ±0.02 


0.99 


± 0.001 


SW Vir 


11.44 ±0.17 


8.64 ±0.16 


1.32 


± 0.030 


RLeo 


17.97 ± 0.02 


14.72 ± 0.03 


1.22 


± 0.002 



The best fit uniform disk radii, and standard deviations, are 
Hsted in Table|3] The last column is the ratio of the L band and 
K band diameters. It is very close to 1 in the case of Betelgeuse, 
but the difference is 20-30% for the other two stars. 

In the next sections, we investigate the possibility that dust 
envelopes around the star may be the cause for the diameter 
variations from K to L. 

3. Modelling of the dust envelope 

The models used in this paper are based on a radiative transfer 
code maintained at the Observatoire de la Cote d' Azur in Nice, 
Fran ce. A descrig tior^^o^^e code is presen ted by Niccolini 
et al. ( l2003l) . We discuss hereafter the main hypotheses and 
parameters of this model. 

3.1. Description of the radiative transfer code 

The code is based on the modelling of stellar radiation and cir- 
cumstellar dust envelope interaction by a Monte Carlo method, 
assuming thermal equilibrium. The star is assumed to be a 
blackbody whose emitted energy is split in a number of pack- 
ets, or "photons". These photons can either pass the envelope 
directly, or they are scattered once or more in the envelope, or 
the photons are absorbed and then thermally re-emitted by the 
dust grains. A grid over the dust envelope is defined with ra- 
dial zones and angular sections, with the resolution accounting 
for expected temperature gradients. From the resulting temper- 
ature profile, the code derives the spectral energy distribution 
(SED) for the three individual photon flux contributions men- 
tioned above. Initially, the stellar radius is left as scaling fac- 
tor. Assigning a specific stellar radius, the total SED can be 
compared to photometric data in order to constrain the model. 
Absorption, emission, and scattering by gaseous components 
are not taken into account. 

Furthermore, the spatial flux distribution at a given wave- 
length can be derived, i.e., a model image of the object at this 
wavelength can be created. The Fourier transform of the model 
image can be compared to visibilities measured at that wave- 
length which yields a further constraint to the model. For the 
purpose of producing the image, the code allows a zooming 
factor while the image dimensions are always fixed at 300x300 
pixels. 

3.2. Parameters of the envelope 

The code offers several parameters that can be chosen. The size 
distribution s(a) of spherical dust particles having the radius a 
was set as s{a) ~ a^^-^, according to the interstellar medium 



jMathis et alJI 19771) and independent of the distance from the 
star. The radius a was contained in the range 0.005 - 0.25 jjm. 

The particle density distribution n{r) of the dust grains that 
was used in this study includes acceleration effects on dust by 
radiation pressure and gas drag (see Appendix IXt. The result- 
ing density distribution is much more peaked at the inner edge 
of the envelope than is a r^^ distribution. 

For our study, the dust particle material was constrained to 
astronomical silicates. Their complex dielectric function (re- 
spectively complex refractive index) is given by Draine & Led 
( 1984). It was sampled at 37 wavelengths between 0.1 jum and 
300 yum, with the sample points considering local extrema. 

The previously mentioned parameters were chosen in con- 
sistency with t he relevant subs et of studies reported bv Danchi 
et al. lil994l) . They were set and fixed in this study. To fit the 
model to the data, the following input parameters were ad- 
justed. The stellar radius and the stellar effective temper- 
ature govern the flux radiated by the star. Spherical symme- 
try was assumed for a dust shell around the star with an inner 
radius n and an outer radius Tq, in units of Furthermore, 
the envelope is characterised by its optical depth, r^, which is 
defined for overall extinction (absorption and scattering). It can 
be chosen at a given wavelength. A, expressed in fim. 

Since the code, as used for this study, does not take into 
account possible changes of the dust properties when heated 
beyond the dust condensation temperature, the temperature 
map has to be checked accordingly for each set of param- 
eters. As condensation temperature for silicates, we adopted 
Tcond ~ 1000 K (Gail&Sedlmavr 1998) . 

4. Comparison between model and observations 

We now study the case of three particular stars. The question 
is whether a simple spherical dust shell model can account for 
the large diameter differences observed in the K and L bands. 
The goal is to find the stellar parameters for which visibility 
points in the K and L bands and at 1 1 fim are compatible and 
simultaneously in agreement with the photometric data. 

4.1. A supergiant: a Orionis 

Betelgeuse is a Semi-Regular (SRc) variable star of spectral 
type MIL It is a late-type supergiant star with a 2335-day pe- 
riod (see Table[3. 

Th e mass-loss rate M of a Orionis is moderate . iKnapp et alJ 
(0^98*) estimate a mass-loss rate from CO emission measure- 
ments of 3.1 • lO^^M© yr ' at HIPPARCOS distance and with 
an outflow velocity Vq = 14.2 km s"'. Throughout this study, 
we adopt terminal outflow velocities from their work. For a 
distance d of 150 p and an outflow velocity Vq = 15 km s"', 
IPanchi et al]lll994h vield a mass-loss rate of 17 ■ lO^^M© yr 
With the HIPPARCOS parallax of 7.63 mas yielding a distance 
of 131 pc, and M oc VqcI, their mass-loss rate is updated to 
14 ■ 10"^ Mq yr"', which is given in Tableg] 

Observational data We used low spectral resolution data 
from different catalogues available at the iSIMBAD. (,2003.) 
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Tabl e 4. Mass-los s rate c omparison. The values for Vq come 
from lKnapp et all ( Il998l) . Data from Danchi et al._ were re- 
evaluated at HIPPARCOS distance (see Table [Hi following 

M ocvod. 



Star 




Mass-loss rate (10 ^ Mq yr ') 




(km s-') 


^Knapp e( al. 


^Danchi e( al. 


^^Present paper 


a Ori 


14.2 


3.1 


14 


17 


SW Vir 


7.8 


1.7 




1.4 


RLeo 


6.8 


0.94 


0.60 


0.77 



Table 5. Parameters adjusted for modelling (see Sect. 13. 2> . T*: 
stellar temperature; 7?*: stellar radius; ro'. inner respectively 
outer radius of dust shell; ti i : optical depth of shell at 11 fim. 
In all cases, an accelerated dust density profile was used. 



Object 




R* 


rJR. 


ro/R^ 


Til 




(K) 


(mas) 








a Ori 


3640 


21.8 


45.9 


48.2 


0.0065 


SW Vir 


2800 


8.65 


15.0 


1000.0 


0.045 


RLeo 


2500 


14.9 


3.5 


140.0 


0.1 



Astronomical Database to constrain the SED over a large range 
of wavelengths. The TDl catalogue provides fluxes around 
0.2 yum. UB V measurements are reported in the GEN and UBV 
catalogues. Fluxes in the UBVRIJHKL bands are taken from 
the JPl 1 Johnson catalogue. The 12, 25, 60, and 100 u m fluxes 
are from the IRAS catalogue, and lDanchi et alJ (1 19941) provide 
flux curves from 8 to 22 fj.m. 

For Betelgeuse, having slight luminosity variations, we can 
easily compare visibility points from different epochs. This as- 
sumption will be fully justified by the fit we obtain. 

Parameters Danchi et al. ('1994') propose a spherical dust shell 
model to account for ISl data in the thermal infrared. This shell 
is supposed to have a small extent (n - 45.97?*, Tq - 48.27?*) 
quite far away from the central star and to be optically very 
thin (tii - 0.0065). It seems to indicate an episodic mass-loss 
process, leading to an empty area between the star and its en- 
velope. We used these parameters (see Table |5jl to reproduce 
the shorter wavelength visibility points given by FLUOR and 
TISIS. 

The effective stellar temperature of 3640 K used bv Danchi 
et al . (,1994.) is fully consistent with the evaluation bv Perrin 
et al. ( l2003al) giving 3640 K and 3690 K depending on models. 
The stellar radius used in that model is 21.8 mas, satisfactorily 
accounting for both the photometric data and interferometric 
data in the three bands. This radius is close to the uniform disk 
model radii in the K and L bands of 21.6 mas and 21.4 mas, 
respectively. 

Results The SED of the model is shown in Fig.^ Apart from 
the total flux from the object, the plot also contains the contri- 
butions of the stellar flux attenuated by the circumstellar shell, 
the distribution of the flux scattered by the dust grains in the 
shell (maximum at short wavelength), and the thermal emis- 
sion of the grains (maximum at longer wavelength). As the 



SED shows, stellar radiation is prominent over a large wave- 
length range. Flux emitted by dust becomes prominent only in 
a small band around 10 jum. 

Visibility curves are also plotted in Fig. ^ The K and L 
band visibility curves from the dust shell model resemble each 
other and show good superposition with the respective visibil- 
ity curves of the uniform disk model. Apparently, the same 
structure is seen in these two bands. The SED plot indicates 
that this should be the stellar photosphere since in these two 
bands the direct stellar flux is roughly 10^ times larger than the 
scattered part and about 10^ times larger than the dust emission. 
One can conclude that the dust envelope, having a very small 
optical depth, has little influence on the intensity distribution 
at these wavelengths. The dust shell is transparent in K and L 
and the star is seen through the dust at these wavelengths. As a 
consequence, the dust shell model derived from the ISl data is 
fully consistent with the FLUOR and TlSlS measurements. 

As a further consistency check, yet with minor priority, we 
computed a mass-loss rate M from our model parameters by 
a method similar to Danchi et al. ( 1994) (see Appendix IaI for 
details). With the set of parameters for a Ori, we obtai n a mass- 
loss rate of 17 ■ 10"^ Mq yr"', consistent with Danchi e^^ 
(^ 1994) and roughly a factor of 5 higher than that of Knapp e t alJ 

Although a Ori has the lowest optical depth of the studied 
objects, its mass-loss rate seems to be the highest. The param- 
eters used here indicate an episodic event in the past. Yet, since 
this event the dust shell has expanded, with its enlarged inner 
radius entering the computation. For a fu rther discussion of the 
mass-loss event see lDanchi et al J ill 9941 

We now study the case of stars with circumstellar envelopes 
at higher optical depths. 

4.2. A semi-regular variable: SW Virginis 

The pulsation period of this M7111 semi-regular (SRb) variable 
star is 150 days. Its main characteristics are recapitul ated in 
Table[2 The mass-loss rate is estimated bv Knapp et al](ll998l) 
to 1.7 ■ lO^^M© yr-i. 

Observational cJata Photometric data are retrieved from the 
IRC catalogue for the 1 and K bands, from the UBV catalogue 
in the respective fi lters, from the IRAS catalogue at 12, 60, 
100 jum, and from iMonnier et al t (.1998) i n the 8 to 13 um 
range. The M band flux at 5 jum comes from lWannier & Sahail 

To date, no 11 -Aim inter f eromet ric data is available for 
SW Virginis. iBenson et aP lll989l) report observations by 
speckle interferometry in the N band at spatial frequencies be- 
low 2-10^ rad Since they do not provide measured data di- 
rectly, but rather give some model parameters, this has not been 
considered here. It is therefore not possible to assess the con- 
sistency of the spatial model between all bands for this object. 

Parameters For SW Vir. lPerrin et alJ ( 1 19981) give an effective 
temperature of 2921 + 100 K, as derived from the bolomet- 
ric flux and the limb-darkened diameter. Consistent with this 
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Fig. 1. Model output for a Ori. From top left to bottom left: modelled visibility curve with interferometric data for the bands K, L, 
and N; bottom right: modelled spectral energy distribution with photometric data. Model parameters (see Table|5ji: = 3640 K, 
= 21.8 mas, n =45.97?*, ro = 48.2/?*, and m =0.0065. 



value, the effective temperature was set to 2800 K. For the pa- 
rameters of the dust shell we c hose to adopt n - 1 5 .0 and 
To — 1000/?*, as proposed bvl van der Veen et alJ (Il995 ^ on 
the basis of sub-mm observations. Applying an optical depth 
Til = 0.045 and assuming a stellar radius of 8.65 mas then 
accounts for flux data at low spectral resolution. 

Danchi et al. ( 1994) make a distinction between two classes 
of objects. The first class has a cold circumstellar envelope far 
from the star (beyond 10 stellar radii), while objects of the sec- 
ond class are surrounded by a warmer dust envelope, near the 
star (at about 2 or 3 stellar radii). The mass-loss process is dif- 
ferent for these two classes. For the second class, the mass-loss 
process is more continuous than for the first one. SW Virginis 



has been classified in the past as a Mira star, hence with a some- 
what regular mass-loss. Yet, its current classification as a semi- 
regular variable, possibly a precursor of a Mira star, makes its 
mass-loss more likely to be irregular, consistently with its ir- 
regular photometric variations. It is therefore very tempting to 
identify it as a member of the first class. The preferred set of pa- 
rameters obtained for this star is summarised in Table|5] A large 
optical depth and a large shell extent are necessary to reproduce 
the photometric contribution of dust at large wavelengths. 



Results As Fig.|2shows, the contribution of dust emission to 
the SED in the K and L bands of SW Virginis is more impor- 
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Spatial Frequency (10^ red"') Wavelength (^tm) 



Fig. 2. Model output for SW Vir. Top left and right: modelled visibility curve with interferometric data for the K and L bands; 
bottom left: expected visibility curve for the N band; bottom right: modelled spectral energy distribution with photometric data. 
Model parameters (see Table||}: = 2800 K, = 8.65 mas, n = 15.0 ro = 1000/?*, and m = 0.045. 



tant than in the case of a Ori, i.e., it is higher in relation to the 
direct flux. This can be attributed to the larger optical depth. 
SW Virginis has some infrared excess that requires a dust shell. 
However, the dust cannot be too warm, otherwise it would pro- 
duce a much larger infrared excess than is compatible with the 
SED. 

From this model, visibility curves were derived. Although 
comparatively few visibility measurements are available for 
this object, it is obvious that the single dust shell model is not 
simultaneously compatible with the K and L band data. Again, 
model visibility curves are similar in K and L which lets us as- 
sume that we see the same structure, i.e. the stellar photosphere. 
This is also reflected by the SED for which a similar reasoning 



applies like in the case of a Ori. Indeed, the dust being quite 
far away from the photosphere, the K band data are well repro- 
duced by the visibility model. However, the dust is too cold and 
optically too thin to lead to a significant increase of the diam- 
eter at larger wavelengths compatible with the L band data. A 
simple dust shell model, therefore, cannot account for both the 
interferometric and the photometric data. Unfortunately, no ISI 
data are available for this star. They would have been very use- 
ful to help model the dust shell and place stronger constraints 
on it. 

The mass-loss rate we obtain from our parameters is 1 .4 • 
10~^ M(T, yr~', which is roughly consistent with the one of 
lKnappetaL(.1998.) . 
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4.3. A Mira-type variable: R Leonis 

R Leonis is a Mira-type long-period variable with a period of 
~ 310 days and of late spectral type MSIIIe ( see Ta ble [p. Its 
mass-loss rate is indicated by Knapp et al. (1998) as 0.94 ■ 
10"^ Mq yr"', whereas the estimate by TPanchiet "ail (11994 
can be updated to 0.60 • lO^^M© yr~'. 

R Leonis is known to have a very strong infrared excess 
and very steady pulsations which makes it very different from 
SW Virginis. R Leonis shows large photometric variations, es- 
pecially at short wavelengths. The problem here is to synthe- 
sise a low spectral resolution SED, as well as visibility curves, 
for a star with large photometric variations. Nevertheless, we 
decided to include all available data points (both photometric 
and interferometric) into the plots despite the object variability 
which will allow a larger range of acceptable model parame- 
ters. 

Observational data For R Leonis, UBV measurements can be 
found in the GEN and UBV catalogues. The JPl 1 Johnson cat- 
alogue provides maanitudes in the BVRIJKL bands. Strecker 
et al. (.1978.) observed this object between 1 .2 and 4.0 fim. The 
fluxes at 12, 20, 60 and 100 um are from IRAS, ari d fluxes be- 
tween 8 and 20 fim come from lDanchi et al.1 (Il994h . 

Parameters The eff'ective temperature that accounts for the 
photometric data is about 2500 K, a bit more than the tem- 
perature of 2000 K assumed bv lPanchi et alJ (11994 *) for their 
model containing only silicate dust around R Leo. The stel- 
lar radius is clearly smaller (14.9 mas instead of 19.8 mas). 
However, we have estimated these two parameters with all data 
combined (interferometric data in the K and N bands and the 
spectro-photometric data), and they show good consistency (in 
as much as the star is neither a perfect blackbody nor a uni- 
form disk) wi th those derived fro m the K band interferometric 
data alone bv lPerrin et alJ (Il999l) . who give uniform disk radii 
of at most ~ 15 mas and photosphere efi'ective temperatures 
closer to 3000 K rat her than 2000 K. Our eff'ective temperature 
is larger than that of lPanchi et al J which is not surprising since 
in their case only photometric data above 8 fim was used, i.e., 
in a range of wavelengths very sensitive to cool dust emission. 
This also explains why their stellar radius is larger: for a given 
bolometric flux, the stellar diameter increases with decreasing 

temperature. 

As in iDanchietalJ (11994 . we consider a large circumstel- 
lar envelope near the central star. But instead of 2 for the 
inner edge of the dust shell, we chose 3.5 iJ* to keep the tem- 
perature of the grains close to 1000 K, the assumed temperature 
of grain condensation. The overall parameters are summarised 
in Tabled 

Results The fit of the photometric data described in the pre- 
vious section is displayed in Fig. |3] The contribution of the 
envelope dominates beyond 10 /im up to about 30 fim. With a 
quite large extent, up to 140 7?*, and an optical depth n i = 0. 1 
the dust sheU accounts for IRAS and IPanchi et alJ (Il994{) 
observations in the long wavelength range. As can be no- 



ticed in the graph, scattered flux and thermal dust emission 
are even more important than for SW Vir. This is apparently 
due to the larger optical depth and the smaller inner radius 
of the envelope. These parameters lead to a mass-loss rate of 
0.77 ■ 10"^ Mq yr"\ which is of the sa me order of and inter- 
me diate between the va lues published bv lDanchietll] ( 11994 
andlKnap p et'ai] ( ll998h (see TablegJ. 

Visibility curves gained from this model are also shown in 
Fig- El for the K, L, and N bands respectively. The stellar ra- 
dius derived from our model is 14.9 mas, much smaller than 
the 19.8 mas of Danchi et al. ( 1994). Since our study is based 
on interferometric data spanning a wider range of wavelengths, 
the parameters we derive are different from that of lPanchi et alJ 
(1994) The ISI data are better suited to characterise the dust, 
yet the shorter-wavelength data bring more accurate stellar pa- 
rameters, hence explaining the discrepancies between the two 
studies. The dust envelope structure may be considered as sta- 
tionary with little sensitivity to stellar pulsation. FLUOR obser- 
vations at phase 0.28 (M arch 1997) in K b and give a uniform 
disk radius of 15.34 mas JPerrin et alJl999l) . In principle, the K 
band measurements are more sensitive to regions very close to 
the photosphere. It therefore seems that the radi us of R Leonis 
should be smaller than the evaluations by Danch iTt alJ (Il994 
and the difference between our estimate and that of Danchi 
et al. cannot be attributed to the star pulsation. This difference 
must be attributed to the difference of eff'ective temperatures. 

We now compare interferometric data in the K and L band 
and that of ISI. Our parameters are still consistent with the ISI 
visibility points at 1 1 yum and account for FLUOR observations 
in the K band at low spatial frequencies. It is to be noticed that, 
since the ISI visibility data are at low spatial frequencies, they 
are not as sensitive to the different models as are the visibilities 
at shorter wa velengths. T his explains why they are consistent 
with both the lPanchi et al.. (.1994.) and our model. On the other 
hand, in the L band our model overestimates the measured vis- 
ibility values. As in the case of SW Vir, it appears that the large 
structure seen in the L band is not well accounted for by a sim- 
ple dust shell model, despite the even larger optical depth. Also 
the SED shows that contributions from the dust shell are at least 
a factor 10^ below the direct stellar flux. A larger radius could 
explain both 3.6 fim and 11 yum observations, but would dis- 
agree with K band measurements and the SED. 

5. Discussion 

5.1. Obtained results 

From what we have seen so far, the best consistency between 
the multiple band data and the dust shell model is obtained for 
a Oh. The model described bv lDanchi et alJ (Il994 has a stellar 
radius close to the K and L radii estimated with a uniform disk 
model. Yet, in this case the influence of the envelope is not 
very important because of its small thickness an d low opacity. 
Nevertheless, model parameters for ISI data from Da nchi et all 
are consistent with what we see both for K and L ba nds. 

One may note, however, that lWeiner et al ] (l2000t) report fur- 
ther ISI data o n a Ori which are not fully consistent with the 
model of iDanchi et al J based on observations at low spatial fre- 
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Fig. 3. Model output for R Leo. From top left to bottom left: modelled visibility curve with interferometric data for the bands K, L, 
and N; bottom right: modelled spectral energy distribution with photometric data. Model parameters (see Table|5ji: 7^* = 2500 K, 
= 14.9 mas, n = 3.5 R*, = 140 R*, andm = 0.1. 



Quencies. lWeiner et alJ took their measurements at higher spa- 
tial frequencies where the dust shell is fully resolved. Its contri- 
bution to the visibility is therefore negligible and the observed 
data are governed by radiation coming directly from the star or 
from its close environment. Li ke the L band d ata of S W Vir and 
R Leo (see below), the data of lWeiner et al .l are overestimated 
by the d ust shell model. These findings are further discussed in 
IPerrin et alJ ( l2003al) . 

For AGB stars with later types like SW Vir and R Leo, for 
which we measure large diameter differences between 2.15 /vm 
and 3.77 yum (the approximate effective wavelengths of the 
FLUOR and TISIS instruments), a simple dust shell model can- 
not account for these differences. However, the shape of near- 



infrared visibility curves for SW Vir and R Leo seem to show a 
trend. L band visibility points are systematically overestimated 
by a simple dust shell model meaning the source is larger than 
predicted by the model. This shows that the stellar environment 
mainly seen at this wavelength is much more complex than sug- 
gested by a model that takes into account only the star and its 
dust envelope. 

Indeed, the temperature of this area is close to the tempera- 
ture at which grains form. Consequently, we must also take the 
presence of gas into account in modelling. Doing this will rad- 
ically change the radiation transfer behaviour. Recently, spec- 
troscopic detection of SiO, OH, and H2O masers in late-type 
star envelopes have revealed the necessity of modelling the 
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gas to explain large diameters seen in L band. ISO/SWS and 
other s tudies (Tsuii et al. ( 1997); Yamamuraet al. ( 1999); Tsuii 
brought up the existence of molecular layers of warm 
and cold H2O, SiO, and CO2 around AGB stars. Such layers 
in AGB stars certainly exist and have an impact on the stellar 
radiation which is observed by infrared interferometers. Strong 
molecular scattering by H2O and CO were suggested to ac- 
count for K band visibihty m easurements of R Leo by FLUOR 
in 1997 (Per rin et alJ ll999> and seem to explain more easily 
the near-infrared measureme nts MMennesson et alJ2002t Perrin 
et al. 2003b). It is very likely that dust needs to be coupled with 
gas to produce a more reliable model of these sources and that 
dust alone cannot be the only cause of the measurements taken 
in a wide range of wavelengths. The inclusion of a detailed pho- 
tospheric model with line opacity due to gas would improve the 
agreement of the SED, particularly in the blue, and would im- 
prove the discrimination between stellar and shell components. 

5.2. Data used 

Certain aspects of our study can be further weighed. 
Simultaneous interferometric data would certainly help con- 
strain models more accurately. Same epoch measurements are 
needed because of the large diameter variations of AGB stars. 
At least for the position angles it can be noted, though, that 
they are fairly close in range, which leaves little possibility that 
the discrepancy we see may be explained by an elliptical shape 
(see below). Furthermore, 1 1 . 15-;um data would have been very 
useful to see whether our SW Vir dust model parameters are 
correct since the envelope has more impact on the model at this 
wavelength. Such data might be provided in the near future by 
single-dish instruments on 10-m class telescopes or by upcom- 
ing mid-infrared interferometers like MIDI at the VLTI. 

5.3. Model code parameters 

A limitation of the model as used for this study is that it does 
not consider heterogeneous grain size distributions s{a) de- 
pending upon the distance from the star or temporal effects 
(like grain formation, growth, or destruction). Also, different 
limits of the grain radius a might have affected the fit of the 
SED in the UV and the far-infrared. Yet, varying the grain size 
range from 0.005-0.25 yum to 0.005-2.5 fim and 0.05-2.5 /im 
for R Leo as a study case, the resulting diagrams showed 
only small changes. It appeared that, in the modelled visibil- 
ity curves for the K and L bands, the contribution of spatial 
frequencies below 20 ■ 10^ rad ' were somewhat increased. 
Above all, neither the star nor the envelope are ideal black- 
bodies which adds to these deviations. 

C oncerning the particle density distribution n(r), ISuhl 
( 1199 9) questions the necessity of assuming a distribution which 
is peaked at the inner edge. Indeed, comparing results of cal- 
culations with two sets of parameters that differ only in n{r) 
showed only minor changes in K and L bands. For 11.15 fim, 
changes were somewhat clearer, whereby for a Ori and R Leo 
the accelerated version fitted ISI data better. Other issues are 
the possible existence of smaller structures in the envelope like 



clumps or aspherical d ust distributio ns, as they were investi- 
gated, for example, by iLopez et al ] (E997). These have not 
been tested here and require more intense visibility modelling. 
This would be definitely required if sources showed system- 
atic asymmetric effects, yet which are likely to produce ringing 
features on top of a visibility curve compared to the more ba- 
sic effect of radiative transfer in spherical shells. Nevertheless, 
closure phase observations could help test for asymmetries. 

As for the dust particle material, slightly differing optical 
properties w ere proposed by several authors. See, for example, 
ISuhl (Il999l) for a comparison of optical properties of silicate 
dust gr ains in the en velopes around AGB stars. We chose the 
ones of lDraine & LeCit 1984) which are w idely used. Using op- 
tical properties suggested bv lSuhl ill999h for comparison made 
it necessary to choose different dust shell parameters in order 
to achieve reasonable consistency with ISI data. Yet, L band 
data were still overestimated, similar to Fig. |3 Another issue 
might be the dust composition. For example, for R Leo, Danchi 
et al. C1994) tested also a mixture of silicates and graphites. 
Yet, this led to a required stellar radius in the dust shell model 
almost twice as large as observed in the K band, which is sup- 
posedly sensitive to the region close to the stellar photosphere. 
Therefore, this has not been further investigated here. A third 
aspect might be the condensation temperature T^ond that we 
adopted. Some authors chose higher Tcond, see, for example, 
Willson (2000) for an overview. Taking R Leo as study case, 
we set ri - 2Ri, and chose and an optical depth ti i so that 
also ISI data were well fitted. The output of the model showed 
a somewhat increased IR excess around 1 1 fim, and the tem- 
perature at the inner radius turned out to be about 1400 K. 
Still, while the SED, K band, and ISI data were reasonably well 
fitted, the modelled visibility curve in the L band was clearly 
higher than the measured data, similar to Fig.|3 

In summary, we do not believe the aspects mentioned above 
to affect our conclusions significantly. 

6. Conclusions 

High angular resolution techniques give important constraints 
on spatial intensity distribution of late-type stars. Using diam- 
eter measurements in the K and L bands, we have tried to ex- 
plain both K and L observations for three stars surrounded by 
circumstellar envelopes with different optical depths. We ap- 
plied a model with a single dust shell and consistent with low 
resolution photometric data as well as thermal infrared inter- 
ferometry data from ISI when available. 

For the supergiant Betelgeuse, the optically very thin dust 
shell has little influence on what is seen and is transparent in 
the K and L bands. The observed diameters are almost equal. In 
contrast, the AGB stars SW Vir and R Leo show bigger diame- 
ters in the L band than in the K band. Here, the optical depths 
of the shells are larger and near-infrared model visibilities react 
quite sensitively on changes of the dust shell parameters. Yet, 
the models we find in agreement with the SED, K band interfer- 
ometry, and the specially dust-sensitive 1 1 yum interferometry 
cannot reproduce the observed diameter differences. 

Our study leads to the conclusion of the insufficiency of a 
simple shell model with dust only or without more complex 
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structures. New observations at several wavelengths fri m vis- 
ible to far-infrared would help understand better external en- 
velopes of AGB stars. A denser coverage of the spatial fre- 
quency plane will also be very useful to strongly constrain 
models. Interferometric observations with higher spectr il reso- 
lution could enable to localise molecular components c f these 
envelopes and certainly help explain differences between K and 
L band. 



Appendix A: Evaluation of the mass-loss rate 

Let Mdust be the constant mass-loss rate of dust ejected from 
the star, ma the average single dust grain mass, and i(r) the 
dust particle density, then at the inner edge ri of the du^t shell 
it holds that 



Mdust = ■ — J- = v'iTOdWi, 

at 



(A.l) 

n(ri) the fcarticle 



where v; is the dust outflow velocity and n 
density. 

In the dust shell, the particle density of grains is assumed 
as fSchutte & Tielens 1989: Danchi et al. 1994.1 



n(r) = m ■ fir) = m ■ ^ 



1 



r2 (l-A-ri/r)i/2 



where the parameter 
A = 1 - (vi/v,f 



(A.2) 



(A.3) 



can be adapted to different cases. It is zero, if Vj is a sumed 
to equal the terminal outflow velocity Vq of the dust far away 
from the star. In this case, the particle density simply vollows 
n(r) oc r^^. For taking into account acceleration of the jiust by 
radiation pressure and coupling with the gas, Vo will be higher 
than Vi. The parameter A will then be larger than in the first 
case and the particle density much more pronounced ^t n. In 
the present study we chose A = 0.99, which applies for typical 
conditions J Panchi et al J 19941) . 

The optical depth of the envelope is determined lly 



f(r)dr. 



(A.4) 



where n(r) was replaced by iA.lt . Cext is the average extinc- 
tion mean section of the dust grains. Replacing ni in €\A\ by 
expression JA.4> yields a relation for the mass-loss rate by dust: 



4nrfvimdTA 
CextW £"f(r)dr' 



(A.5) 



In the previous two expressions, r\, Tq, and are the cho- 
sen model parameters. Cext is calculated in the modelling code 
for a given grain size distributi on from the op tical properties 
of the dust by Mie theory ( Niccolini et al. "2003V Likewise, the 
integral in ( IA.5I I is performed by the code. The dust grain mass 
ffJd has to be estimated by averaging the size distributi|)n s(a) 
and defining the mass density. For the presently considered sil- 
icate dust, we chose as pd = 3.3 g cm"^ (also used bv iDraine 



& Leel ( 1 1984 and lOanchi et alJ ( Il994 ). The terminal veloc- 
ity of he shell particles i s taken from high- spectral resolution 
observations of CO lines ('Knapp etalJl998h . 

Furthermore, for the total mass-loss rate also the gas has to 
be taken into account, i.e.. 



Md„„ + M„; 



(A.6) 



Gas-tc-dust mass ratios are not well known and differ quite a 
lot. W; adopted here, as typical value and assuming dust and 
gas to have the same terminal velocities jSchutte & Tielen^ 
Il989: Danchi et aL1994:.Habing.l996) . 



250 Mdu 



(A.7) 



This means that most of the mass-loss is carried by gas. For 
Mgas, an expression analogous to ( lA.U can be given where Vo 
enters for Vi. Considering iA.ll . and therefore neglecting the 
mass-l'oss carried by dust, the total mass-loss rate is eventually 



M, 



total 



250 



AnrhonidTA 
CM J'°fir)dr' 



(A.8) 
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